Synthesis, Characterization, Fluorescence Properties, and DFT Modeling of Difluoroboron Biindolediketonates

We report a simple and efficient strategy to enhance the fluorescence of biocompatible biindole diketonates (bdks) in the visible spectrum through difluoroboronation (BF2bdks complexes). Emission spectroscopy testifies an increase in the fluorescence quantum yields from a few percent to as much as >0.7. This massive increment is essentially independent of substitutions at the indole (-H, -Cl, and -OCH3) and corresponds to a significant stabilization of the excited state with respect to non-radiative decay mechanisms: the non-radiative decay rates are reduced by as much as an order of magnitude, from 109 s−1 to 108 s−1, upon difluoroboronation. The stabilization of the excited state is large enough to enable sizeable 1O2 photosensitized production. Different time-dependent (TD) density functional theory (DFT) methods were assessed in their ability to model the electronic properties of the compounds, with TD-B3LYP-D3 providing the most accurate excitation energies. The calculations associate the first active optical transition in both the bdks and BF2bdks electronic spectra to the S0 → S1 transition, corresponding to a shift in the electronic density from the indoles to the oxygens or the O-BF2-O unit, respectively.


Introduction
Difluoroboron β-diketonates (BF 2 bdks) are a class of tetracoordinate organoboron complexes that has been investigated only in the last 16 years, after the first example reported by Fraser and co-workers in 2007 [1]. These complexes have shown excellent photophysical properties, such as broad UV-Vis absorption range, large extinction coefficients, two-photon excited fluorescence, fluorescence in solution up to the NIR region, and, in some cases, solid-state luminescent behavior and room-temperature phosphorescence (RTP).
Another captivating feature of these adducts is the ability to produce singlet oxygen ( 1 O 2 ), which allows BF 2 bdks to be used as oxygen sensors in tissue visualization or as photosensitizers (PSs) in photodynamic therapy (PTD) [2].
Some modifications in molecular conformation and intermolecular interactions in the solid state influence the luminescence of BF 2 bdks, making them responsive to external stimuli, such as pressure, heat, or mechanical force. Thermochromic or mechanochromic fluorescence is useful in sensoristics for high-end applications [3,4].
Substituents at C1 and C3 of the diketone ligand greatly influence the behavior of the boron compounds, and slight structural modifications on the diketone ligands enable the fine tuning of luminescence properties. Extended π conjugation is known to positively affect fluorescence emission efficiency [5], and intramolecular charge transfer (ICT) is especially favored when the electron-poor dioxadifluoroborine ring portion is coupled with electron-rich aryl scaffolds, along with methoxy or amino groups [6].
Despite their extended π conjugation, bdks are not strongly luminescent for two main reasons: Firstly, the structures of the diketo tautomers are non-planar. Secondly, in the keto-enol tautomer, the activation of very fast, non-radiative decay mechanisms is provided by intra-and intermolecular interactions of the acidic enol proton [18]. Here, it is shown that the bdk fluorescence quantum yield can be greatly enhanced by reacting bdks with boron trifluoride etherate: this enables it to lock the 1,3-diketone into its enolate form and to introduce an acidic, electron-poor center onto an electron-rich backbone. Time-dependent (TD) density functional theory (DFT) calculations have been used to investigate the electronic structure of these new BF2bdks, with the aim to assign the peaks observed in the experimental spectra and to explain (i) the red shift of the maximum in the experimental UV-Vis absorption spectra upon difluoroboronation, and (ii) the increase in the fluorescence lifetime of BF2bdks with respect to the corresponding bdks. Additionally, we show that these new BF2bdks combine the fluorescence properties of the BF2bdk core and its ability to generate singlet oxygen with the indole moiety, which, notably, is bioactive.

Synthesis
Friedel-Crafts acylation, a type of electrophilic aromatic substitution, is a wellknown and viable strategy for the functionalization of aromatic and heteroaromatic Scheme 1. Reaction pathway.
Despite their extended π conjugation, bdks are not strongly luminescent for two main reasons: Firstly, the structures of the diketo tautomers are non-planar. Secondly, in the keto-enol tautomer, the activation of very fast, non-radiative decay mechanisms is provided by intra-and intermolecular interactions of the acidic enol proton [18]. Here, it is shown that the bdk fluorescence quantum yield can be greatly enhanced by reacting bdks with boron trifluoride etherate: this enables it to lock the 1,3-diketone into its enolate form and to introduce an acidic, electron-poor center onto an electron-rich backbone. Timedependent (TD) density functional theory (DFT) calculations have been used to investigate the electronic structure of these new BF 2 bdks, with the aim to assign the peaks observed in the experimental spectra and to explain (i) the red shift of the maximum in the experimental UV-Vis absorption spectra upon difluoroboronation, and (ii) the increase in the fluorescence lifetime of BF 2 bdks with respect to the corresponding bdks. Additionally, we show that these new BF 2 bdks combine the fluorescence properties of the BF 2 bdk core and its ability to generate singlet oxygen with the indole moiety, which, notably, is bioactive.

Synthesis
Friedel-Crafts acylation, a type of electrophilic aromatic substitution, is a well-known and viable strategy for the functionalization of aromatic and heteroaromatic substrates. It allows one to obtain aryl ketones in a single synthetic step starting from a suitable aryl partner and an acyl chloride, catalyzed via a Lewis-acid-like aluminum trichloride. If the aromatic partner is electron-rich, addition of the acid is optional.
With the idea to combine the indole moiety with the 1,3-diketone nucleus (see Scheme 1), we selected three different indole substrates as starting materials (1 in Scheme 1): 5H-indole (a), 5-chloroindole (b), and 5-methoxyindole (c). These substrates were chosen as they have a neutral (a), a slightly electron-poor (b), and an electron-rich molecule (c), respectively. The three commercial NH indoles were alkylated at the nitrogen atom (step 1 → 2 in Scheme 1) to protect it from undesired electrophilic attack and unwanted boron complexation, in addition to conferring solubility to end complexes. The use of methyl iodide and sodium hydride as the base in dry THF under an inert atmosphere gave the methylated products in high yields (intermediates 2 in Scheme 1). Diketone synthesis (intermediates 3) was accomplished using malonyl dichloride, without the addition of AlCl 3 , as already reported by Cirrincione and co-workers [14]. Good yields were obtained for 5H-N-methylindole (3a) and 5-methoxy-Nmethylindole (3c), while poor results were obtained with 5-chloro-N-methylindole (3b). To verify if the poor results obtained with the chlorinated compound were associated to the electron poverty of the substrate, we repeated the acylation using an electron-poorer substrate, 5-nitro-N-methylindole. In this case, no conversion of the starting material at all was observed. Accordingly, in the literature, there are very few examples of acylation of electron-poor indoles, always requiring very strong Lewis acids such as stannic chloride. The synthesized diketones are labeled in the following ways: HBIP (compound 3a), HClBIP (compound 3b), and HBMIP (compound 3c).
The isolated diketones were then treated with excess boron trifluoride etherate in dichloromethane (step 3 → 4), causing the precipitation of the corresponding organoboron complex in the reaction mixture after a few minutes of stirring at room temperature. Filtration and washing of the precipitates usually gave the products in a satisfactory yield and good purity for the following spectroscopic analyses. The BF 2 bdks complexes are labeled BF 2 BIP (compound 4a), BF 2 ClBIP (compound 4b), and BF 2 BMIP (compound 4c).
The as-synthesized BF 2 bdks 4 are bright-yellow powders, with an evident fluorescence under a 365 nm UV lamp, both in the solid state and in solution. Compounds bdks 3 and BF 2 bdks 4 were characterized using 1 H and 13 C-NMR spectroscopy as well as infrared (IR) spectroscopy (see Supplementary Materials). Further investigation was carried out with 11 B and 19 F-{ 1 H}NMR spectroscopy (see Supplementary Materials), confirming the structure and purity of the synthesized BF 2 bdks.
Considering the keto-enol tautomerism, each bdk examined in this study is a mixture of three conformers: keto-enol, trans-diketo, and cis-diketo (see Figure A1c, Figure A1b, and Figure A1a, respectively). The relative amount of these conformers is influenced by different factors, mostly the temperature and the medium (solvents of different polarity, bulk versus gas phase . . . ). 1 H-NMR indicates a larger concentration of diketo with respect to keto-enol in all the solvents, although an increase in the latter is observed as the solvent polarity rises (see Table 1). −60 cm −1 upon difluoroboronation (∆ ν C=O , see Figure S1), indicating coordination at the carbonyl group (see Table S1). These results are confirmed via the DFT calculations, predicting a redshift going from 3a to 4a (see Section S5 and Figure S2).

UV-Vis Absorption Spectroscopy
The UV-Vis absorption spectra of the BF 2 bdks 4a, 4b, and 4c were recorded in a panel of solvents differing as to their polarity and hydrogen bonding properties, namely: toluene (aprotic, non-polar solvent), ethyl acetate, dichloromethane, acetone, acetonitrile (polar solvents, weakly H-bonding acceptors), dimethylformamide, dimethyl sulfoxide (strong H-bond acceptors), butanol, ethanol, and methanol (protic solvents). For details on the sample preparation and measurement apparatuses, see Section 3.4. To evaluate the effects of fluoroboronation on the spectroscopic properties of the compounds, similar measurements were also performed on the corresponding non-modified compounds 3a, 3b, and 3c (see Figure S3).
In Figure 1, we report the UV-Vis absorption spectra of 4a (a), 4b (b), and 4c (c) in toluene, dichloromethane, dimethyl sulfoxide, and methanol. The wavelengths corresponding to the maxima of the absorption spectra (λ abs ) in all the solvents considered are reported in the third column of Table 2. The value of the molar extinction coefficient at λ abs , ε, as determined through linear fit of absorbance vs. concentration plots, is also reported in the fourth column of Table 2 (see Section 3.4 for details). Table 2. UV-Vis absorption properties of 4a, 4b, and 4c (see Figure 1). Maximum absorption wavelength (λ abs ) and molar extinction coefficient at λ abs (ε).
Like other enolized β-diketones [19][20][21][22][23][24][25][26], the BF 2 bdk compounds show, in all the solvents, a main absorption band peak in the blue portion of the visible spectrum, between 428 nm (4b in toluene) and 452 nm (4c in dimethyl sulfoxide), and a much fainter absorption in the UVB (around 280 nm). The spectral line shape of the main peak exhibits a substructure, with a shoulder around 400 nm. A systematic redshift is observed for all three compounds on increasing the solvent polarity, and a further bathochromic effect is induced by H-bonding interactions with H-bond-accepting solvents. Nevertheless, λ abs is only slightly dependent on both the environment and the phenolic substituents (see Table 2), unlike what is observed for other classes of BF 2 complexes [27]. The molar extinction coefficient ε varies by more than a factor of 2 on both the solvent and the phenolic substituent, although a systematic dependence cannot be evidenced. However, the absorption cross-section is notably high in any environment for all the compounds.
In terms of variance, the UV-Vis absorption spectra of the non-fluoroboronated compounds 3a, 3b, and 3c (see Figure S3) are characterized by three bands, two in the UV region, peaking around 255 and 305 nm, and another one in the violet, between 350 nm and 400 nm (see Table 3). Table 3. UV-Vis absorption properties of 3a, 3b, and 3c. Peak absorption wavelengths in each of the three absorption bands resolved within the electronic transition spectra: λ abs1,diketo , λ abs2,diketo , and λ abs,enol . f diketo : fractional concentrations of trans-diketo (see text for details on calculation) 1 . The maximum absorption wavelengths measured for each of the non-fluoroboronated compounds (3a, 3b, and 3c) in all the solvents are detailed in Table 3. Based on the calculations (see Section 2.3.1), we attribute the UV peaks to absorption of the trans-diketo conformer (λ abs1,diketo and λ abs2,diketo in Table 3) and the remaining peak (λ abs,enol in Table 3) to the enol conformer. The blueshifted peak of the enol species suggests reduced conjuga-tion in the non-modified compounds with respect to the fluoroboronated ones. The typical redshift from UV to the blue region of λ abs upon difluoroboronation [10] is, thus, also confirmed for the present class of compounds. The presence of relevant amounts of compounds in the cis-diketo conformers, which was hypothesized for other β-diketones [23,24], is excluded in DFT calculations in this instance (see Appendix A).

Compound
By assuming that the value ε(λ abs2,diketo )/ε(λ abs,enol ) = 0.356, calculated using TD-B3LYP-D3 for 3a in toluene (see Section 2.3.1), does not vary appreciably upon either changing the solvent or shifting from 3a to the other two compounds, we can estimate the fraction of compound in trans-diketo conformation as: where Abs(diketo) and Abs(enol) represent the absorbances detected at λ abs2,diketo and λ abs,enol , respectively, in all solvents but acetone. The trans-diketo fractional concentration values, calculated according to Equation (1), are reported in Table 3. This analysis confirm that observed with 1 H NMR (see Table 1), that the keto-enolic equilibrium is shifted towards the diketo species in all the cases in which it can be estimated (i.e., in any solvent but acetone).

Steady-State Fluorescence Measurements
The fluorescence emission spectra of the fluoroboronated compounds 4a, 4b, and 4c, were recorded in all the solvents upon excitation at the pertaining λ abs . The emission spectral line shapes in toluene, dichloromethane, dimethyl sulfoxide, and methanol are plotted in Figure Table 4. In the same table, the fluorescence quantum yield values, φ fluo , determined via comparison with dimethyl-POPOP dissolved in cyclohexane (φ fluo = 0.95, [26]), are also reported (see Section 3.4 for details on the calculation of φ fluo ).
The emission appears to be poorly solvent-dependent, the only appreciable solvation effect being a sizeable redshift from the non-polar solvent toluene to any of the others.
With the only exception of 4c in methanol (φ fluo = 0.07), fluorescence is extremely intense in any of the tested solvents. Namely, the fluorescence quantum yield value is >0.17 in all the environments.
The compounds 3a, 3b, and 3c emit fluorescence in distinct bands upon excitation at either of the absorption peaks. For all three species, the fluorescence of the diketo conformer is negligibly low. This might be due to a drastically reduced conjugation extent because of the non-planarity of the molecular structure. Excitation at λ abs,diketo resulted in comparably faint emission peaking around 340 nm. The fluorescence quantum yield values could barely be estimated and were lower than 0.01 in all of the tested solvents. The fluorescence peak wavelengths and quantum yields for 3a, 3b, and 3c measured upon excitation at the enol absorption peak, λ abs,enol , are listed in Table 5. The emission of the enol fraction of the unmodified compounds is systematically blueshifted with respect to that of the corresponding fluoroboronated species in all the solvents. Moreover, the fluorescence quantum yields, although being notably higher with respect to those measured for the diketo tautomers, which reflects enhanced conjugation in the planar enol conformers, are typically more than one order of magnitude lower than those measured in the case of the fluoroboronated analogues.     The excitation spectra of the fluoro-boronated compounds were also recorded, with the observation wavelength set at 500 nm. Exemplary spectra are shown in Figure 3. The spectral line shapes, although somewhat redshifted, are very similar to those of the corresponding UV-Vis absorption spectra, suggesting very simple excited-state dynamics.

Time-Resolved Fluorescence Measurements
In order to investigate the excited-state dynamics of both the non-substituted compounds 3 and their fluoroboronated analogues 4 in depth, time-resolved fluorescence measurements were performed by exploiting a time-correlated single-photon counting apparatus endowed with 30 ps temporal resolution. For the fluoroboronated species, the decay is single exponential in all the tested solvents. An exemplary decay pattern is plotted in Figure 4. The best-fitting curve is also shown as a solid line. The fluorescence lifetimes of BF 2 bdks in the different solvents are reported in Table 4. In agreement with steady-state fluorescence data, the values of the fluorescence lifetime are high and rather independent from variations in both the polarity and the H-bonding properties of the solvent. Moreover, 4a, 4b, and 4c behave very similarly. This pattern suggests very simple decay dynamics, dominated by radiative decay and basic photophysical mechanisms, such as internal conversion and intersystem crossing. Combining the lifetime data with the measured quantum yield values and exploiting Equations (2), we could estimate the radiative (k Fl ) and non-radiative (k NR ) decay rates, whose values are reproduced in Table 4.   in the activation of efficient non-radiative decay mechanisms, with the overall effect of incrementing kNR by roughly two orders of magnitude.
According to the above-proposed rationale, the dramatic increase in fluorescence quantum yield observed upon fluoroboronation may be ascribed to the concomitant hindrance of photoinduced tautomerization and excited-state intramolecular proton transfer dynamics.

Singlet Oxygen Production
The generation of 1 O2 upon irradiation with UV-Vis light is the main source of the photobiological activity of a class of drugs named Type II photosensitizers. Upon electronic excitation, these pharmaceutical active principles undergo intersystem crossing to the metastable T1 state lying beneath the fluorescent singlet and the electronic ground state. As the molecular oxygen ground state is a spin triplet, collisional quenching between a photosensitizer and oxygen results in a selection-rules-allowed concomitant transition of the drug to the ground state and of oxygen to its spin-singlet radical form. This latter is extremely reactive. Accordingly, if produced within biological tissues, it is capable of inducing severe oxidative stress, thereby triggering apoptosis as well as other cell death mechanisms. Because the ability of a putative drug substance to act as Type II photosensitizer ultimately depends on the relative probability to decay through intersystem crossing, stabilizing the S1 state with respect to other non-radiative decay pathways may be beneficial to the optimization of the photopharmaceutical performance of a molecule. In this view, as a preliminary evaluation of the potential of the compounds of class 4 as photosensitizers, the 1 O2 generation quantum yields, ϕ( 1 O2), of solutions of 3a, 3b, 3c, 4a, 4b, and 4c were assessed relative to that of meso-tetraphenylporphyrin tetrasulphonate (TSPP), a well-characterized photosensitizer currently undergoing clinical trials [28]. The amount of produced 1 O2 was estimated in terms of the oxidation-induced quenching of the 9,10-dimethylanthracene fluorescence [29]. Singlet oxygen quantum yields are reported in Table 6 (see Section 3.6 for details on the experimental procedures).
While none of the bdk compounds exhibit significant 1 O2 generation efficacy, the corresponding BF2bdks all display sizeable photosensitizing activity, as expected based on what was reported previously for other BF2bdks [2]. Although the activity of compounds The above calculations suggest that the high luminescence of the compounds is due to the concomitance of a fast radiative decay rate with a comparable (thus, unusually slow) non-radiative decay rate. Moreover, also these indirect parameters are rather independent from both environmental properties and substitutions at the lateral aromatic rings.
The situation is much more complicated for compounds of class 3. Fluorescence decay time distributions were reconstructed both upon excitation at 280 nm (in the absorption spectral region typical of the diketo tautomer) and upon excitation at 420 nm (in the absorption band of the enol conformers). Details on the excitation sources and timing setup are given in Section 3.4. Upon excitation in the UV, biexponential decays were invariably obtained. The values of decay times, τ i , and relative amplitudes, f i , as retrieved from fitting of the experimental decay patterns (see Section 3.4 for details on the fitting procedure), are reported in Table S2. A long-lived component with time constant ∼ = 2 ns, similar to that detected for the fluoroboronated compounds, is still observed, but the dominant transient is a much shorter one, in the range 200-400 ps, which we attribute to a non-radiative decay mechanism encompassing photoactivated transition to the enol structures, in analogy to what was observed in other β-diketones [23,24]. The situation is very similar upon excitation at 420 nm (See Table S3), where the 200-400 ps component, ascribable to the well-established reketonization process, occurs with even higher relative amplitude with respect to long-lived radiative decay (i.e., has relatively higher probability to take place than transition from the diketo to the enol structures). This reflects the superior thermodynamic stability of diketo with respect to enol conformers, indicated by UV-Vis absorption data (see Table 3), 1 H NMR (see Table 1), and demonstrated via DFT calculations (vide Appendix A). Moreover, for all the compounds of class 3, a very short decay time, not exceeding a few tens of picoseconds, is detected in several solvents. We interpret this transient as due to non-radiative decay through exchange of the enol proton between the keto oxygens. The fact that this mechanism is not always evidenced might relate to a relatively weak ketoenolic intramolecular H-bonding affinity, which might, in turn, explain the unusual diketo-shifted tautomeric equilibrium of these compounds. Calculation of the average excited-state lifetime, τ Av , according to Equation (3), allows one to apply Equation (2) to estimate k Rad and k NR .
The estimated rate constants indicate that, although k Rad is somewhat reduced in the bdks with respect to BF 2 bdks, the main reason is that the fluorescence quenching reseeds in the activation of efficient non-radiative decay mechanisms, with the overall effect of incrementing k NR by roughly two orders of magnitude.
According to the above-proposed rationale, the dramatic increase in fluorescence quantum yield observed upon fluoroboronation may be ascribed to the concomitant hindrance of photoinduced tautomerization and excited-state intramolecular proton transfer dynamics.

Singlet Oxygen Production
The generation of 1 O 2 upon irradiation with UV-Vis light is the main source of the photobiological activity of a class of drugs named Type II photosensitizers. Upon electronic excitation, these pharmaceutical active principles undergo intersystem crossing to the metastable T 1 state lying beneath the fluorescent singlet and the electronic ground state. As the molecular oxygen ground state is a spin triplet, collisional quenching between a photosensitizer and oxygen results in a selection-rules-allowed concomitant transition of the drug to the ground state and of oxygen to its spin-singlet radical form. This latter is extremely reactive. Accordingly, if produced within biological tissues, it is capable of inducing severe oxidative stress, thereby triggering apoptosis as well as other cell death mechanisms. Because the ability of a putative drug substance to act as Type II photosensitizer ultimately depends on the relative probability to decay through intersystem crossing, stabilizing the S 1 state with respect to other non-radiative decay pathways may be beneficial to the optimization of the photopharmaceutical performance of a molecule. In this view, as a preliminary evaluation of the potential of the compounds of class 4 as photosensitizers, the 1 O 2 generation quantum yields, φ( 1 O 2 ), of solutions of 3a, 3b, 3c, 4a, 4b, and 4c were assessed relative to that of meso-tetraphenylporphyrin tetrasulphonate (TSPP), a well-characterized photosensitizer currently undergoing clinical trials [28]. The amount of produced 1 O 2 was estimated in terms of the oxidation-induced quenching of the 9,10-dimethylanthracene fluorescence [29]. Singlet oxygen quantum yields are reported in Table 6 (see Section 3.6 for details on the experimental procedures). Table 6. Singlet oxygen quantum yields, φ( 1 O 2 ), of the bdk and BF 2 bdk compounds dissolved in ethanol relative to TSPP(aq).
0.066 * From Ref. [28]. ** The fluorescence of 9,10-dimethylanthracene is quenched less than in the absence of putative photosensitizers after irradiation with the Wood lamp.
While none of the bdk compounds exhibit significant 1 O 2 generation efficacy, the corresponding BF 2 bdks all display sizeable photosensitizing activity, as expected based on what was reported previously for other BF 2 bdks [2]. Although the activity of compounds 4 is lower than that of TSPP, these results confirm the effectiveness of the fluoroboronation strategy to achieve excited-state stabilization and, thereby, the enhancement of intersystem crossing dynamics.

DFT Calculations
Cartesian coordinates for all optimized geometries can be found as XYZ files in the Supplementary Materials.

UV-Vis Absorption Spectra
All the functionals adopted in the study provided a qualitatively identical description of the modification of the UV-Vis spectra upon difluoroboronation and functional group on the indole. The best quantitative agreement between the computed and experimental spectra was obtained using B3LYP-D3 without TDA. For this reason, the results obtained with this level of calculation are discussed later, while the results obtained with the other methods are reported in the Supplementary Materials (see Figures S4-S8 and Tables S4-S8). This result was unexpected based on a recent assessment of several density functionals for the prediction of excitation energies of most organic and main-group molecules [30]. In this assessment, ωB97X-D was the functional characterized by a lower root mean square error in the prediction of the excitation energies, while B3LYP performed poorly. Here, the trend is the opposite. Interestingly, although the HOMO-LUMO band gap (∆E HOMO-LUMO ) is a rough approximation of the excitation energy, its value computed with B3LYP-D3 almost coincides with the excitation energy (see Table 7), while all the other functionals largely overestimate it (see Tables S4-S8). This result agrees with what was reported for other compounds [31,32]. For what concerns the TDA approximation, it introduces a negligible error with respect to TD, underestimating λ abs of an additional 10-15 nm. Because of the lower computational cost of TDA calculations, it is confirmed to be a suitable approximation for the modeling of excitation energies also for bdks and their complexes. Table 7. Excited state properties of 3a keto-enol, 3a trans-diketo, 4a, 4b, and 4c. Vertical (∆E vert ), adiabatic (∆E adia ), and ZPVE-corrected adiabatic (∆E adiaZPVE ) excitation energies for the S 0 → S 1 transition in the gas phase. The HOMO to LUMO energy difference in S 0 (∆E HOMO-LUMO ) and the vertical de-excitation energy for S 1 → S 0 transition (∆E vert ) are also showed for comparison. All the energies are reported in eV. The absorption wavelength (λ absB3 , nm) corresponding to the S 0 → S 1 transition is also showed with its oscillator strength (f absB3 , a.u.), ground to excited state transition electric dipole moment ( S0→S1 a.u.), and the radiative decay rates of the S 1 state (k radB3 , s −1 ), defined as in Equation (4). k radB3 is computed considering the intensity and the maximum of the corresponding band in the computed spectrum. The emission wavelength (λ fluoB3 ) corresponding to S 1 → S 0 computed on S 1 geometry, and its oscillator strength (f fluoB3 ) is also reported. All the values have been obtained at DFT/TD-DFT with B3LYP-D3/def2-TZVPD. The computed UV-Vis absorption spectra of the three conformers of 3a are reported in Figure 5a, while the spectra of trans-diketo 3a, 4a, 4b, and 4c are shown in Figure 5b. For what concerns the three conformers of 3a, their computed spectra are significantly different for wavelengths > 270 nm: the keto-enol has a maximum at 362 nm (see Table 7), while the trans-diketo and the cis-diketo peak at 309 and 283 nm, respectively. Based on the relative energy stability of the three 3a conformers, both in the ground and first excited state (see Appendix A, Tables A1 and A2), we expect to observe signals associated only to trans-diketo and keto-enol conformers in the experimental spectra. The computed spectra for the keto-enol and the trans-diketo match the signals recorded experimentally in toluene (see Table 3 and Figure S3), supporting the assignments discussed in Section 2.2.1. While the band at 362 nm (oscillator strength fabsB3 = 0.70) of the keto-enol is determined to be associated to a single transition, the band at 309 nm of trans-diketo is the convolution of three signals at 331.34 (fabsB3 = 0.08), and two almost degenerate at 314.30 (fabsB3 = 0.02) and 314.28 nm (fabsB3 = 0.14). The signal observed in the experiments at ~245 nm is associated based on the calculations to a higher-energy electronic transition of transdiketo (excited state 12, 247.9 nm, fabsB3 = 0.23).
The difference between the total electronic density of the first excited state and of the ground state is shown for 3a trans-diketo (331 nm) and 3a keto-enol (362 nm) in Figure  6a,b, respectively. As the oscillator strength of the transition to the third excited state is double that of the first excited state for trans-diketo, the corresponding total density difference is also reported in Figure 6a'. For both conformers, the electronic transitions correspond to a transfer of the charge from the indole and the CH2/CH group of the 1,3dicarbonyl to the C=O/C-O(H). While for 3a trans-diketo, the transition causes a charge transfer at short range, often on the same atom (see, for example, the oxygen atoms in Figure 6a,a'), for 3a keto-enol, the charge transfer involves different atoms, i.e., in the excited state, the regions with a lower electron density than the ground state (orange background in Figure 6b) are localized on different atoms than the regions with a higher electron density (blue background). The larger charge displacement is at the origin of the larger ground to excited state transition electric dipole moment (  ) of the first electronic transition in 3a keto-enol than 3a diketo (see Tables 2 and S4-S8). For what concerns the three conformers of 3a, their computed spectra are significantly different for wavelengths > 270 nm: the keto-enol has a maximum at 362 nm (see Table 7), while the trans-diketo and the cis-diketo peak at 309 and 283 nm, respectively. Based on the relative energy stability of the three 3a conformers, both in the ground and first excited state (see Appendix A, Tables A1 and A2), we expect to observe signals associated only to trans-diketo and keto-enol conformers in the experimental spectra. The computed spectra for the keto-enol and the trans-diketo match the signals recorded experimentally in toluene (see Table 3 and Figure S3), supporting the assignments discussed in Section 2.2.1. While the band at 362 nm (oscillator strength f absB3 = 0.70) of the keto-enol is determined to be associated to a single transition, the band at 309 nm of trans-diketo is the convolution of three signals at 331.34 (f absB3 = 0.08), and two almost degenerate at 314.30 (f absB3 = 0.02) and 314.28 nm (f absB3 = 0.14). The signal observed in the experiments at~245 nm is associated based on the calculations to a higher-energy electronic transition of trans-diketo (excited state 12, 247.9 nm, f absB3 = 0.23).
The difference between the total electronic density of the first excited state and of the ground state is shown for 3a trans-diketo (331 nm) and 3a keto-enol (362 nm) in Figure 6a,b, respectively. As the oscillator strength of the transition to the third excited state is double that of the first excited state for trans-diketo, the corresponding total density difference is also reported in Figure 6a'. For both conformers, the electronic transitions correspond to a transfer of the charge from the indole and the CH 2 /CH group of the 1,3-dicarbonyl to the C=O/C-O(H). While for 3a trans-diketo, the transition causes a charge transfer at short range, often on the same atom (see, for example, the oxygen atoms in Figure 6a,a'), for 3a keto-enol, the charge transfer involves different atoms, i.e., in the excited state, the regions with a lower electron density than the ground state (orange background in Figure 6b) are localized on different atoms than the regions with a higher electron density (blue background). The larger charge displacement is at the origin of the larger ground to excited state transition electric dipole moment (µ S0→S1 ) of the first electronic transition in 3a keto-enol than 3a diketo (see Table 2 and Tables S4-S8). After functionalization with BF2, the experimental spectra show the typical redshift of the absorption maxima from the UV to the visible region, as reported in the literature [10]. The TD-DFT spectra reproduce the shift in the absorption maximum, going from 3a (keto-enol, dark-violet curve in Figure 5b) to 4a (yellow curve). The computed spectra for compounds 4 show the same line shape of the 3a keto-enol spectrum, although redshifted (see Figure 5b). The maximum absorption wavelength in toluene is experimentally observed at 431 nm for 4a, while it is slightly blueshifted to 428 in 4b and slightly redshifted to 436 nm for 4c (see Table 2). These changes are qualitatively supported by the calculations, with corresponding computed values of 386.6 (4a, R = −H, fabsB3 = 0.88), 386.2 (4b, R = −Cl, fabsB3 = 0.94), and 405.2 nm (4c, R = −OCH3, fabsB3 = 0.85), respectively. This small effect of the R group on the excitation energies is supported by the CM5 charges computed on the BF2 unit and the 1,3-dicarbonyl that are very similar in all the BF2bdks. The experimental band observed at ~430 nm for compounds 4 is associated based on the TD-DFT calculations to a single transition. The shoulder at 400 nm is assigned alternatively to a vibronic profile or to local distortion of the orientation of the two indoles from the full planarity (see Section 2.3.2). Figure 6c shows, for 4a, the effect on the electronic total density of the electronic excitation corresponding to the experimental 431 nm signal. This plot is essentially identical to the one obtained for 3a, showing a charge transfer from the indole to the carbonyls and the fluorine. The computed oscillator strengths for compounds 4 are larger than those computed for the corresponding bdk compounds. While for the keto-enol forms, the increase in fabsB3 is only 25%, for the trans-diketo, i.e., the majority species in the experiments, it is 500%. The larger fabsB3 in compounds 4 can be explained by the larger charge displacement occurring in the excitation, going, in order, from trans-diketo (Figure 6a,a') to ketoenol (Figure 6b) to BF2 complexes (Figure 6c). Based on Figure 6, a larger displacement of After functionalization with BF 2 , the experimental spectra show the typical redshift of the absorption maxima from the UV to the visible region, as reported in the literature [10]. The TD-DFT spectra reproduce the shift in the absorption maximum, going from 3a (ketoenol, dark-violet curve in Figure 5b) to 4a (yellow curve). The computed spectra for compounds 4 show the same line shape of the 3a keto-enol spectrum, although redshifted (see Figure 5b). The maximum absorption wavelength in toluene is experimentally observed at 431 nm for 4a, while it is slightly blueshifted to 428 in 4b and slightly redshifted to 436 nm for 4c (see Table 2). These changes are qualitatively supported by the calculations, with corresponding computed values of 386.6 (4a, R = −H, f absB3 = 0.88), 386.2 (4b, R = −Cl, f absB3 = 0.94), and 405.2 nm (4c, R = −OCH 3 , f absB3 = 0.85), respectively. This small effect of the R group on the excitation energies is supported by the CM5 charges computed on the BF 2 unit and the 1,3-dicarbonyl that are very similar in all the BF 2 bdks. The experimental band observed at~430 nm for compounds 4 is associated based on the TD-DFT calculations to a single transition. The shoulder at 400 nm is assigned alternatively to a vibronic profile or to local distortion of the orientation of the two indoles from the full planarity (see Section 2.3.2). Figure 6c shows, for 4a, the effect on the electronic total density of the electronic excitation corresponding to the experimental 431 nm signal. This plot is essentially identical to the one obtained for 3a, showing a charge transfer from the indole to the carbonyls and the fluorine. The computed oscillator strengths for compounds 4 are larger than those computed for the corresponding bdk compounds. While for the keto-enol forms, the increase in f absB3 is only 25%, for the trans-diketo, i.e., the majority species in the experiments, it is 500%. The larger f absB3 in compounds 4 can be explained by the larger charge displacement occurring in the excitation, going, in order, from trans-diketo (Figure 6a,a') to keto-enol ( Figure 6b) to BF 2 complexes (Figure 6c). Based on Figure 6, a larger displacement of the charge is obtained in BF 2 bdks than in bdks, even considering the keto-enol conformers because of the involvement of the fluorine atoms. This results in larger µ S0→S1 and then oscillator strengths: µ S0→S1 computed for 3a (keto-enol) and 4a are, in fact, 9.4 and 11.2 a.u., respectively (see Table 7).

UV-Vis Emission Spectra
The optimized geometry for the first excited state of all the compounds is very similar to the one obtained for the ground state at M06-2X, ωB97X-D, and CAM-B3LYP-D3. TDA-B3LYP-D3 and TD-B3LYP-D3 instead predicted that, as result of the first excitation, the compound structure goes from a geometry having the two indole coplanars in S 0 (as in Figure A1c) to a geometry very similar to trans-diketo in S 1 for 3a keto-enol, 4a, 4b, and 4c using TDA-B3LYP-D3 and 4c using TD-B3LYP-D3 (see Figure S11). The coplanar conformers in S 1 are computed to be less stable by~20 kJ mol −1 , while the other levels of calculations agree in predicting the two indoles coplanar also in S 1 . For what concerns 3a cis-diketo, the first excited state for these functionals will undergo cyclization. These results are independent on the number of states included (up to 100), the maximum of the optimization step, and the presence/absence of D3 correction. The experimental emission and excitation spectra support the similarity of the ground and excited state geometry for all the compounds and then a planar geometry for 3a keto-enol, 4a, 4b, and 4c in S 1 . The computed emission spectra for a non-planar geometry for compounds 4 would be significantly different than the corresponding absorption spectra (see the emission spectra of the two 4c conformers in Figure 7). Reoptimization of the non-planar conformers of 3a keto-enol, 4a, 4b, and 4c in S 1 using M06-2X, ωB97X-D, and CAM-B3LYP-D3 restores the planar geometry. Based on this assessment, B3LYP-D3 seems to not be a suitable method for the optimization of excited states for the bdks and BF 2 bdks class. Nevertheless, TD-B3LYP-D3 level provides the transition energies closest to the experimental values, as already noted for the absorption spectra. In order to facilitate a comparison with the experiments, the results reported in the following for S 1 were obtained using TD-B3LYP-D3, using singlepoint calculations for 3a cis-diketo and for 4c on the  Electronic fluorescence spectra of (a) the three conformers of 3a and (b) the BF2bdks compounds as computed at the TD-B3LYP-D3/def2-TZVPD level including the first 100 excitations. Two spectra are reported for 4c: the spectrum obtained on the S1 geometry optimized with TD-B3LYP-D3 (4c no planar, see Figure S11) and with TDA-CAM-B3LYP-D3 (4c planar).
All the methods agree with the experiments in the qualitative description of the effect of difluoroboronation and of the chemical modification of the indole on the emission spectra. The calculations quantitatively reproduce the shift in the emission bands with respect to the absorption spectra for all the compounds. While for 3a, a redshift is predicted of 53-16 cm −1 , for the BF2 complexes, the shift is about +20 cm −1 , independently of the R group. The wavelengths corresponding to the maximum of the emission spectra in the >350 nm region are also very similar for all the BF2bdks, as verified for the absorption spectra: 406.0 Two spectra are reported for 4c: the spectrum obtained on the S 1 geometry optimized with TD-B3LYP-D3 (4c no planar, see Figure S11) and with TDA-CAM-B3LYP-D3 (4c planar).
All the methods agree with the experiments in the qualitative description of the effect of difluoroboronation and of the chemical modification of the indole on the emission spectra. The calculations quantitatively reproduce the shift in the emission bands with respect to the absorption spectra for all the compounds. While for 3a, a redshift is predicted of 53-16 cm −1 , for the BF 2 complexes, the shift is about +20 cm −1 , independently of the R group. The wavelengths corresponding to the maximum of the emission spectra in the >350 nm region are also very similar for all the BF 2 bdks, as verified for the absorption spectra: 406.0 (4a, R = −H, f absB3 = 0.88), 407.1 (4b, R = −Cl, f absB3 = 0.94), and 422.9 nm (4c, R = −OCH 3 , f absB3 = 0.87), respectively (see also Figure 7 and Figures S4-S8, Table 7 and Tables S4-S8). Again, very similar charges are computed on the BF 2 unit and the 1,3-dicarbonyl also in the excited state, indicating a poor influence of the R group on the electronic properties of the compounds. The computed radiative rate constants (k rad , see Table 7 for B3LYP-D3 and Tables S4-S8 for the other methods) are able to reproduce the increase of one order of magnitude in the experimental k Rad upon difluoroboronation (see k Rad in Tables 4 and 5). Diketo and the keto-enol conformers of bdks have very different computed electronic properties, with the keto-enol having an extinction coefficient and radiative decay rate associated with the first excitation energy one order of magnitude larger than diketo ones. Thus, the increase in the experimental k Rad upon difluoroboronation is a further confirmation of the larger concentration of diketo with respect to keto-enol conformers in the solutions. In fact, k radDFT for 3a keto-enol is only slightly lower than the 4a one, and both are one order of magnitude larger than 3a trans-diketo. This result also explains the increase in k Rad with the solvent polarity in compounds 3 (see Table 4), as associated with the increase in the relative concentration of keto-enol, as actually quantified with 1 H NMR (see Table 1).
The increase in the fluorescence properties of BF 2 bdks can then mainly be associated with the stabilization of the keto-enol geometry. A positive contribution is also brought about by the introduction of the BF 2 unit (although less important than the stabilization of a planar and more symmetric geometry), and it is associated with a decrease in the band gap with respect to the keto-enol conformers and an increase in the transition electric dipole moment. The increase in the dipole moment is explained by total density maps, showing a larger charge displacement in the first electronic transition in BF 2 bdks than in bdks because of the involvement of the BF 2 unit (see Figure 6). Confirming the larger role played by the symmetry of the geometry in enhancing the fluorescence properties of these compounds, the lowest-energy band of the emission spectrum of a distorted 4c conformer (see Figure  S11 Table 7).

Materials
All reactions were performed in oven-dried glassware under normal atmosphere, unless otherwise noted. All reagents and solvents were purchased from Fluorochem, TCI, or Merck and used without any further purification. Where specified, reactions were monitored via thin-layer chromatography (TLC) on POLYGRAM ® Xtra Sil G/UV254 (0.2 mm layer thickness; Macherey-Nagel); spots were observed under a UV lamp at 254 or 365 nm. Gravimetric column chromatography was performed using silica gel (60 Å, particles size: 0.63-0.2 mm) as stationary phase. 1 H, 13c , 19 F, and 11 B NMR spectra were acquired on a Bruker AVANCE I 400 instrument at 400.13, 100.61, 376.5, and 128.38 MHz, respectively, and are referenced using residual non-deuterated solvents (CHCl 3 7.26 ppm in CDCl 3 , DMSO 2.50 ppm in d 6 -DMSO); deuterated solvents were purchased from Eurisotop; multiplicities are abbreviated as singlet (s), doublet (d), doublet of doublets (dd), and triplet (t); coupling constants are reported in Hz; IR spectra were acquired in attenuated total reflectance (ATR) mode using a FT-IR Thermo Scientific Nicolet iS10 Smart iTR instrument equipped with a diamond optical element, over a range 4000-650 cm −1 (at
Synthesis of N-methyl indoles 2b and 2c and diketones HBIP 3a, HClBIP 3b, and HBMIP 3c followed a reported literature procedure (see Ref. [14] and Section S1 in the Supplementary Materials). tiplicities are abbreviated as singlet (s), doublet (d), doublet of doublets (dd), and triplet (t); coupling constants are reported in Hz; IR spectra were acquired in attenuated total reflectance (ATR) mode using a FT-IR Thermo Scientific Nicolet iS10 Smart iTR instrument equipped with a diamond optical element, over a range 4000-650 cm −1 (at 4 cm −1 resolution). Intensities are denoted as: br = broadened signal, sh = shoulder, vs = very strong, s = strong, m = medium, w = weak and vw = very weak.

Materials
All reactions were performed in oven-dried glassware under less otherwise noted. All reagents and solvents were purchased fr Merck and used without any further purification. Where specifie tored via thin-layer chromatography (TLC) on POLYGRAM ® Xt layer thickness; Macherey-Nagel); spots were observed under a U Gravimetric column chromatography was performed using silica 0.63-0.2 mm) as stationary phase. 1 H, 13c , 19 F, and 11 B NMR spe Bruker AVANCE I 400 instrument at 400.13, 100.61, 376.5, and 12 and are referenced using residual non-deuterated solvents (CH DMSO 2.50 ppm in d6-DMSO); deuterated solvents were purchas tiplicities are abbreviated as singlet (s), doublet (d), doublet of d (t); coupling constants are reported in Hz; IR spectra were acqu reflectance (ATR) mode using a FT-IR Thermo Scientific Nicolet iS equipped with a diamond optical element, over a range 4000-65 tion). Intensities are denoted as: br = broadened signal, sh = shou = strong, m = medium, w = weak and vw = very weak.

1,3-Bis(1-methyl-3-indolyl)propane-1,3-dione Difluorobora
To a solution of HBIP (100 mg, 0.30 mmol) in dichloromethan mL, 9.0 mmol) is added, and the mixture is stirred at rt for 5 m forms. The solid is filtered on a Buchner filter funnel and was methane and then dried under vacuum. BIP difluoroborate is obt (73 mg, 64% yield). 1  about by the introduction of the BF2 unit (although less important than the stabilizat a planar and more symmetric geometry), and it is associated with a decrease in the gap with respect to the keto-enol conformers and an increase in the transition elect pole moment. The increase in the dipole moment is explained by total density showing a larger charge displacement in the first electronic transition in BF2bdks th bdks because of the involvement of the BF2 unit (see Figure 6). Confirming the large played by the symmetry of the geometry in enhancing the fluorescence properties of compounds, the lowest-energy band of the emission spectrum of a distorted 4c conf (see Figure S11) is associated with five excitations (λfluoB3 in nm, ffluoB3 in a.u.): (674.95, (450.25, 0.03), (409.13, 0.31), (366.22, 0.01), and (349.96, 0.32). All these bands have a ffluoB3 than the corresponding planar conformer (0.87 a.u., see Table 7).

Materials
All reactions were performed in oven-dried glassware under normal atmospher less otherwise noted. All reagents and solvents were purchased from Fluorochem, T Merck and used without any further purification. Where specified, reactions were tored via thin-layer chromatography (TLC) on POLYGRAM ® Xtra Sil G/UV254 (0. layer thickness; Macherey-Nagel); spots were observed under a UV lamp at 254 or 36 Gravimetric column chromatography was performed using silica gel (60 Å, particle 0.63-0.2 mm) as stationary phase. 1 H, 13c , 19 F, and 11 B NMR spectra were acquired Bruker AVANCE I 400 instrument at 400.13, 100.61, 376.5, and 128.38 MHz, respec and are referenced using residual non-deuterated solvents (CHCl3 7.26 ppm in C DMSO 2.50 ppm in d6-DMSO); deuterated solvents were purchased from Eurisotop tiplicities are abbreviated as singlet (s), doublet (d), doublet of doublets (dd), and (t); coupling constants are reported in Hz; IR spectra were acquired in attenuated reflectance (ATR) mode using a FT-IR Thermo Scientific Nicolet iS10 Smart iTR instru equipped with a diamond optical element, over a range 4000-650 cm −1 (at 4 cm −1 r tion). Intensities are denoted as: br = broadened signal, sh = shoulder, vs = very str = strong, m = medium, w = weak and vw = very weak.

Synthesis of the Complexes
To a solution of HBIP (100 mg, 0.30 mmol) in dichloromethane (100 mL), BF3·Et2 mL, 9.0 mmol) is added, and the mixture is stirred at rt for 5 min: a colored preci forms. The solid is filtered on a Buchner filter funnel and washed with cold dic methane and then dried under vacuum. BIP difluoroborate is obtained as an orange (73 mg, 64% yield). 1

Steady-State Electronic-State Transition Spectroscopy Measurements
The UV-Vis absorption spectra were recorded with a Perkin Elmer Lambda2 spectrophotometer. The fluorescence emission and excitation spectra were recorded with a PTI fluorescence master system spectrofluorimeter. The instrument was interfaced with the acquisition software Felix 2000 (version 1), which performed online correction of the data with respect to the excitation lamp spectral radiance and detector spectral quantum efficiency.
The solvents used to prepare the solutions for all electronic-state transition spectroscopy experiments were purchased from Merck and were of HPLC grade. They were used without further purification. Due to the poor solubility of the compounds in some of the solvents, a concentrated stock in acetone was prepared weighing some powder and dissolving it in 3 mL of this solvent. Afterwards, 50 µL aliquots of this stock were put in glass vials and let to evaporate for 48 h at 50 • C. Finally, the as-obtained films were resuspended in 3 mL of the desired solvent, and the samples were sonicated for 1 h to favor disaggregation. From this procedure, clear solutions were obtained in all solvents.
Molar extinction coefficients were estimated by performing linear regressions on absorbance vs. concentration plots, obtained through progressive dilution of the above stock solutions to one-tenth of the initial concentration in several steps (N > 5). An exemplary plot is reported in Figure S10.
Fluorescence quantum yields, φ Fluo , were determined by comparing the integrated fluorescence intensity with that of dimethyl-POPOP, a dye used in particle physics as a reactant of scintillator chambers. Due to its high-fluorescence quantum yield, which is tabulated to be φ Fluo = 0.95 when dissolved in cyclohexane and excited in its UVA absorption band [26], dimethyl-POPOP has been also widely used as a fluorescence standard since the 1970s. Fluorescence was elicited exciting the samples at their absorption maxima; the emission spectra were normalized with respect to their peak absorbances. Fluorescence quantum yield values were corrected for the solvent's refractive indexes.

Time-Resolved Fluorescence Measurements
The time-resolved fluorescence decay patterns of the compounds in solution were reconstructed referring to the time-correlated single-photon counting (TCSPC) technique. The experimental apparatus exploited to this aim is fully described elsewhere [23,24]. Its overall time resolution can be estimated in terms of the temporal point-spread function (TPSF), i.e., the temporal profile of the excitation laser pulses directly sent to the detectorsensitive area, as reconstructed by the TCSPC system. Typical TPSF at different wavelengths is reproduced in Figure S9. For the experiments reported hereby, the excitation source was a Ti/saffire cw-mode-locked laser (mod. Tiger-ps, Time Bandwidth Products, Zurich, CH), delivering pulses at 48 MHz repetition rate of 3.9 ps duration at the fundamental wavelength. The samples were excited at 420 nm, in proximity of the enol tautomers absorption peak, by the built-in second harmonic of the beam. To elicit the diketo tautomer fluorescence in the non-fluoroboronated compounds 3a, 3b, and 3c, the third harmonic was generated out of cavity, as described elsewhere [23,24]. The fluorescence photons were detected by means of an MPD50 single-photon avalanche diode with integrated active quenching and cooling circuitry and timed by a single module of an SPC-152 integrated board (Becker & Hickl GmbH, Berlin, Germany). The decays were fitted to single-, double-, or triple-exponential model functions exploiting the Levenberg-Marquardt algorithm built in the software Origin 7 for parameter optimization. The number of components of each decay was established adding transients one by one until addition of further decay components did not lead to improvement in the fit quality (evaluated in terms of X 2 value and fit residuals). The data were acquired in triplicate, and the fitting parameters reported in this article are averages over the three parallels. The errors are expressed in terms of the pertaining standard deviations.

Estimation of Singlet Oxygen Generation Efficiency
The improved efficiency in photosensitizing singlet oxygen ( 1 O 2 ) generation of the BF 2 bdk with respect to the corresponding bdk compounds was assessed in ethanol (the most biocompatible solvent in our panel) by exploiting the DMA fluorescence assay [29]. Due to the marked tendency of this sensor to oxidate, a fresh 1 mM concentrated stock of DMA in ethanol was prepared from the powdered compound immediately before each singlet oxygen generation measurement according to the procedure detailed elsewhere [25]. Solutions of the BF 2 bdk at approximate concentration of 1 µM were prepared and their exact concentration determined spectrophotometrically using the molar extinction coefficients reported in Table 2. Due to the impossibility of deriving the precise concentration of bdks in the enol conformers, we chose to assess the 1 O 2 generation of solutions of these latter compounds having the same absorbance value ( ∼ =0.05) at the enol absorption peak. Finally, the ability of BF 2 bdks to quench the DMA fluorescence (i.e., to generate 1 O 2 ) was compared to that of a 1 µM concentrated solution of TSPP in water. The exact concentration was assessed spectrophotometrically assuming the molar extinction coefficient value ε(413 nm) = 129,000 M −1 cm −1 at the Soret band peak [28]. DMA was added from the 1 mM stocks at a final concentration of 20 µM. The obtained samples were placed in 1 × 1 cm 2 fluorimetry quartz cells to carry out the fluorescence measurements. The DMA fluorescence was excited in correspondence of its 340 nm absorption peak and measured in the band 390-650 nm using the PTI spectrofluorimeter described in Section 3.4, prior to and after photosensitization of 1 O 2 through irradiation of the solutions for 30 s with a Wood lamp. Control measurements were carried out in the same experimental conditions on DMA alone.
Geometry optimization was carried out by means of the Berny optimization algorithm with analytical gradient and default convergence thresholds. Unscaled, harmonic vibrational frequencies were computed analytically. Zero-point vibrational energies (ZPVE) were calculated at 1 atm and 298 K from conventional ideal gas, rigid rotator, particle in a box, and quantum mechanical harmonic oscillator partition functions. Charge and spin densities were obtained using Charge Model 5 (CM5) [43] and Hirshfeld population analysis [44], respectively.
Time-dependent (TD) density functional theory in the frame of the linear response theory [45] is the most common approach for the study of molecular electronic excited states [46,47]. All the calculations were performed using the Tamm-Dancoff approximation (TDA, [30,48]). Liang et al. [30] showed that TDA offers a good balance between computational efficiency and accuracy. For comparison, B3LYP-D3 and ωB97X-D calculations have also been performed without TDA (TD-B3LYP-D3). The electronic absorption (emission) spectra were computed by single-point calculations on the ground (first excited)-stateoptimized structures by considering the first 100 excitations. All the calculations were performed on the gas phase single molecule, other than for the evaluation of the relative stability of 3a conformers, where a polarizable continuum toluene model was also adopted (ε = 2.3741 [49]). A Gaussian broadening was applied to the TDDFT excitations using a stan-dard deviation of 0.3 eV. No scaling factors were adopted. Visualization of the computed spectra and of the total density maps was performed using GaussView 6.0.
The vertical excitation energies (∆E vert ) and purely electronic adiabatic excitation energies (∆E adia ) for S 0 → S 1 were obtained as the energy difference between the ground state energy and the energy of the first excited state computed through single-point and geometry optimization, respectively. No difference was observed between geometry optimization, including the first 5 and 100 excitations in an initial assessment using ωB97X-D. For this reason, only the first five excitations were considered in the geometry optimization. ZPVE-corrected vertical adiabatic excitation energies were also obtained (∆E adiaZPVE ). A graphical definition of these energies is reported in Figure 8. The corresponding values for the S 1 → S 0 transition are indicated with a prime (see Figure 8). The radiative decay rates of the S1 species can be estimated as the Einstein coefficient for the spontaneous emission, as defined in Refs. [50,51]: where , c, and ħ represent the vacuum permittivity, the speed of light, and the reduced Planck constant, respectively, while   is the ground to excited state transition electric dipole moment.

Conclusions
In this article, we report on the synthesis of three biocompatible biindole diketonates and their BF2 compounds. Difluoroboronation is confirmed to be an effective strategy to increase, by as much as one order of magnitude, the fluorescence quantum yield of organic compounds, presenting a suitable chelating group [10,27,52,53] and, in particular, for diketonates [6,10]. The above features make these BF2bdks interesting for low-energy-consumption lighting devices, such as organic light-emitting diodes (OLEDs [54]) and fluorescent sensors. Chemical modification of the indole does not significantly influence the electronic spectroscopy of these BF2bdks. The radiative decay rates k rad of the S 1 species can be estimated as the Einstein coefficient for the spontaneous emission, as defined in Refs. [50,51]: where ε 0 , c, andh represent the vacuum permittivity, the speed of light, and the reduced Planck constant, respectively, while µ S0→S1 is the ground to excited state transition electric dipole moment.

Conclusions
In this article, we report on the synthesis of three biocompatible biindole diketonates and their BF 2 compounds. Difluoroboronation is confirmed to be an effective strategy to increase, by as much as one order of magnitude, the fluorescence quantum yield of organic compounds, presenting a suitable chelating group [10,27,52,53] and, in particular, for diketonates [6,10]. The above features make these BF 2 bdks interesting for low-energyconsumption lighting devices, such as organic light-emitting diodes (OLEDs [54]) and fluorescent sensors. Chemical modification of the indole does not significantly influence the electronic spectroscopy of these BF 2 bdks.
These experimental observations were qualitatively reproduced by all the (TD)DFT functionals used. Although TD-B3LYP-D3 showed the best agreement with the experimental excitation energies among the levels used (overestimation of~50 nm), this method was determined to be unreliable for the determination of the excited-state geometries, for which other methods (e.g., ωB97X-D) should be preferred.
The new BF 2 bdks presented in this article combine their biocompatibility with fluorescence properties and the indole moiety, which is notably bioactive. We show that the stabilization of the S 1 excited state towards deactivation in BF 2 bdks is large enough to enable sizeable 1 O 2 photosensitized production. The above figures of merit promise future perspectives of these compounds for applications in biology and medicinal chemistry.

Conflicts of Interest:
The authors declare no conflict of interest.
Sample Availability: Samples of the compounds are available from the authors upon request.

Appendix A
The DFT-optimized structure of the three conformers of HBIP (3a) is reported in Figure A1: trans-diketo, cis-diketo, and keto-enol. Several input structures were used for the three conformers, for which optimization resulted in one of the three conformers reported in Figure A1. Only the keto-enol compound is planar, while the cisand transdiketo prefers a bent geometry, with an angle of~110 • between the rings. This difference is expected based on the different hybridization (sp 2 and sp 3 , respectively) of the C atom connecting the C=O groups in the keto-enol and the diketo conformers, respectively. For the cis-diketo, the two C=O groups are rotated to maximize the distance of the two negatively charged oxygen atoms, i.e., to minimize their repulsive interaction.
The relative energy stability of the three conformers of 3a in the gas phase is listed in Table A1, as obtained using the different functionals here employed. All the levels of approximation provide coincident results and agree in indicating the trans-diketo as the most stable conformer. Its energy is lower by about 10 kJ mol −1 than the keto-enol and 25 kJ mol −1 than cis-diketo. The relative stability of the conformers is slightly changed if the optimization is conducted by including toluene as implicit solvent. Based on a Boltzmann distribution, we can estimate at 25 • C a relative ratio among trans-diketo:keto-enol: cisdiketo of 94:6:0.
We also compared the relative energy stability of the 3a conformers in the first excited state in Table A2. All the levels employed agree that 3a cis-diketo is the most unstable conformer, as already discussed for the ground state. For what concerns the relative stability of 3a trans-diketo and 3a keto-enol, this value is both functional and TD/TDA-dependent, not allowing one to draw clear conclusions like for the ground state. three conformers, for which optimization resulted in one of the three conformers reported in Figure A1. Only the keto-enol compound is planar, while the cis-and trans-diketo prefers a bent geometry, with an angle of ~110° between the rings. This difference is expected based on the different hybridization (sp 2 and sp 3 , respectively) of the C atom connecting the C=O groups in the keto-enol and the diketo conformers, respectively. For the cisdiketo, the two C=O groups are rotated to maximize the distance of the two negatively charged oxygen atoms, i.e., to minimize their repulsive interaction. The relative energy stability of the three conformers of 3a in the gas phase is listed in Table A1, as obtained using the different functionals here employed. All the levels of approximation provide coincident results and agree in indicating the trans-diketo as the most stable conformer. Its energy is lower by about 10 kJ mol −1 than the keto-enol and 25 kJ mol −1 than cis-diketo. The relative stability of the conformers is slightly changed if the optimization is conducted by including toluene as implicit solvent. Based on a Boltzmann  Table A2. Relative energy stability (in kJ mol −1 ) of the first excited state of the three conformers of 3a in the gas phase, as obtained using different TDDFT functionals.